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Abstract. A model based on simple assumptions about 4-dimensional space-time 
being closed and isotropic, and embedded in a 5 th large-scale dimension r representing 
the radius of curvature of space-time, has been used in an application of Newton's 
Second Law to describe a system with angular momentum. It has been found that 
the equations of MOND used to explain the rotation curves of galaxies appear as a 
limit within this derivation and that there is a universal acceleration parameter, a Q , 
with a value, again consistent with that used by MOND. This approach does not 
require modification of Newtonian dynamics, only its extension into a fifth large-scale 
dimension. The transition from the classical Newtonian dynamics to the MOND regime 
emerges naturally and without the introduction of arbitrary fitting functions, if this 
5-dimensional model is adopted. The paper also includes the derivation of an effect 
in 5-dimensional orbital dynamics which is in reasonable agreement with the observed 
Pioneer Anomaly. 



PACS numbers: 04.50+h 



Orbital dynamics in 5 dimensions 



2 



1. Introduction 

A discrepancy exists between the extent of observed visible matter in the universe and 
that inferred from the motion of galaxies and clusters of galaxies if Newton's Second 
Law and Law of Gravitation are applied. The existence of abundant unobserved (or 
"dark" ) matter has been proposed by many cosmologists to account for the discrepancy. 
An alternative but no less radical explanation is that either: Newton's Second Law 
should be modified for the low accelerations evident in observations of galactic motion; 
or that the equations of gravity should be modified for these low accelerations [1,2]. 
MOND was first put forward in 1983 [1] and during the last twenty three years has 
continued to be consistent with the growing body of observational data on the rotation 
curves of galaxies [3], and has survived a number of studies looking for consistency 
with various other effects otherwise requiring dark matter, such as cosmological models 
derived from combining a multitude of observational data [4], to tidal streams from a 
companion galaxy to the Milky Way [5] . However, MOND Theory is empirical in nature, 
being derived exclusively from correlations based upon observations and, to date, has 
eluded derivation from first principles, although the recent relativistic more fundamental 
formulation from Bekenstein [5] offers both MOND and Newtonian limits. 

In this paper we show that the introduction of a 5 th large-scale dimension to the 
description of space-time can account for all aspects of the empirical equations of MOND, 
provide a natural transition function [i between the MOND and Newtonian regimes, and 
explain the link between the MOND parameter a Q and key universal parameters, such 
as Hubble's Constant. An essential feature of this approach is that it does not require 
the principles of Newton's Second Law to be modified, only that they should be applied 
in five dimensions of space-time rather than four. 

The paper is divided into four sections as follows: formulation of Newton's Laws in 
a 5-dimensional space-time; the origin and role of the acceleration, a r , originating from 
the presence of a 5 th dimension; connection to the equations of MOND; and a derivation 
and calculation of the Pioneer Anomaly. The paper concludes with a short discussion 
section. 

2. Formulate Newton's Laws in a 5-dimensional space-time 

The 5 th dimension represented by the parameter r is assumed orthogonal to the three 
space dimensions s(x,y,z) and to the time dimension t. It is also assumed that bodies 
can only propagate in x, y, and z and have no access to the r dimension. However forces 
originating from the r dimension can act on bodies in x, y, and z. This is illustrated in 
figure 1, where only a single space dimension x is shown together with the 5 th dimension 
r forming an x — r plane. 

A body at point o is subject to three forces (accelerations): (i) gravitational 
acceleration g x towards a gravitational source mass M in the direction of negative x; (ii) 
an applied inertial acceleration a x in the direction of x positive; and (iii) an assumed 
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Figure 1. Force balance on a body at point o in the x — r plane. 



(MOND-like) acceleration a r in the positive r direction and orthogonal to x. In a world 
of three space dimensions and no 5 th dimension, balance would simply be defined in 
terms of component balances in each of the three space dimensions s(x,y,z). With a 
5 th dimension it is useful to introduce the concept of branes and look for a condition 
in which the components in two mutually orthogonal branes (denoted by symbols o 
and <7j_) can be balanced. In figure 1 the angle at which the branes intersect with the 
reference axes (x, r) is denoted by 0. 
Force balance occurs when 

in brane a: — g x cos (f> — a r sin + a x cos (1) 
in brane a±: — g x sin + a r cos <fi = a x sin (2) 

These two equations are simultaneously satisfied if 

g x = -a x ± ia r (3) 

Ordinarily it would be assumed that the imaginary part could not be realised by 
measurement in 'real' space. However, as will be shown later this imaginary part can 
give rise to measurable effects under certain circumstances. 

Before leaving this section we return to the consider the angle 0, which defines 
where the brane walls are found. This angle is given by 

tan0 = — (4) 
ax 

This angle will not be constant throughout space but will be determined by the local 
value of this gradient. 
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Figure 2. Relationship between curvature of space-time r u and the acceleration a r . 



3. The origin and role of acceleration, a r , originating from the presence of a 
5 th dimension 

3.1. The origin of a r 

In a closed universe, which is isotropic on a large-scale, there is a finite radius of curvature 
to space-time even remote from gravitating bodies. We attribute this underlying 
parameter, r u , to our assumed 5 th dimension. Propagating electro-magnetic radiation 
effectively follows a world-line of radius of curvature r u and this can be interpreted as a 
universal acceleration applying throughout space (see section 4.1) with a value 



acting in the direction of positive rjj 

The validity of this relationship implies that a closed 4-dimensional space-time 
is considered to be embedded within a locally flat 5-dimensional universe. This is 
illustrated in figure 2 from which we now derive a r . 

In figure 2 the space dimension s represents any one of the three space dimensions 
s(x,y,z). A single space dimension segment Sis 2 represents a distance As. c is 
the speed of light in vacuo, and the radius of curvature of a (geodesic) path of 
electromagnetic radiation in space-time remote from gravitational fields is r u . The 
radius of curvature of single space dimension segment, As, (remote from gravitational 
fields) is then given by r s = r u sin<x Assuming Newton's second law applies in 5 
dimensions, the centrifugal acceleration acting on electromagnetic radiation travelling 
along space dimension segment As is given by: a s = c 2 /r s The value of a r is now given 
by the component of the centrifugal acceleration acting in the r direction - i.e. in the 

X This acceleration has, of course, already been empirically associated with the MOND acceleration 
parameter a Ql through the equation a Q — c 2 / R Cl where R c is the curvature [7] 




(5) 



Orbital dynamics in 5 dimensions 



5 




Figure 3. Branc (<r) defined by balance between gravitational (g x ) and curvature 
accelerations (a r ). 



direction of expansion of the Universe. Hence: 

c 2 c 2 

a r = — sin a = — (6) 



3.2. Considering the role of a r 

At any point o in space dimension x adjacent to a gravitating mass M, but otherwise 
empty, there exist two accelerations; (i) the acceleration due to the gravitational field of 
M, and (ii) a r acting in direction of r increasing. This is shown in figure 3. At the point 
o, a brane can be defined in which there exists balance between these two accelerations, 
g x , and a r . This brane a satisfies the equation: 

fl , a , = -a r tan0 (7) 

where tan0 is the local gradient, dr/dx, of the brane at point o relative to the space 
dimension x. 

The brane (defined as a in the diagram) is the locus of all points o at which there 
is balance between g x and a r , i.e. between accelerations of contraction of the Universe 
on the one hand and expansion of the Universe on the other. It can be convenient to 
think of this balance condition as a form of "surface tension" somewhat akin to the 
"rubber- sheet" analogy of General Relativity. This balance condition is a function of 
distance x, away from the mass, and this allows a local functional form for tan0(x) in 
the presence of a gravitating body as: 

tan 0= ^ = 1^™ = ^: ( 8 ) 

a,,x 2 c 2 x 2 dx 



r (x) = r\ 



GM 



C 2 X 



(9) 



An interesting characteristic of the brane a can be noted at this point simply from 
the equation of balance between g x and a r : the closer is point o to the gravitating mass 




Figure 4. Rate of change in direction of velocity. 



M, the greater is the gradient dr/dx and the more closely aligned is the brane a to 
the r-dimension, and as we move further away and out towards infinity the more the 
acceleration lying in the plane of the brane a itself converges to g x . 



4. Connection to the equations of MOND 

The original motivation for MOND was to explain the quasi-flat rotation curves of 
galaxies. Hence to establish if there is a connection between this work and MOND we 
need to see how the force balance in 5 dimensions affects orbiting systems. For an object 
in a stable circular orbit (i.e. constant velocity, v) in the x — y plane (see figure 4) the 
centrifugal acceleration in standard 4-dimensional space-time is 

a x = Ql z A V (10) 

Adapting this for 5 dimensions requires the introduction of a new angle, 9, between the 
angular momentum vector, Q zr , which now lies in the z — r plane and the r axis as 
shown in figure 5. Equation 10 then becomes 

a>x = 2r A VL — \^zr\ (sin# + i cosd) v (11) 
where the imaginary component is in the r direction. 

The angle 9 will again be locally determined and by contrast with the variation 
in the angle <fi which arises from the dependence of potential energy on position, we 
associate 9 with the kinetic energy. Replacing potential energy with the non-relativistic 
kinetic energy in equation 9 giveq§| 



r (v) = 7\ 



v 2 ' 
0.5^ 

c z 



(12) 



\VL ZT \ is independent of 9 and can be determined from the familiar case of Q zr being 
aligned strictly in the direction of z (i.e. 9 = 7r/2) when |0^ r | = v /(radius of orbit). 
As for 9 we attribute this to the gradient, dr/dx, induced by the inertial accelerations 
acting on a body, and this must be calculated from equation 12 using 

dr dr dv 
dx dv dx 



It is interesting to note that equation 12 suggests we can also write r (v) = r u yf\ — v 2 /c 2 
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Figure 5. Angular velocity vector £l zr lying in the plane zr. 




A r (x+Ax) 



Figure 6. : The orbital plane lies in xyr space at an angle 9 to the xy plane. 



For the body in orbit in figure 4 this gives 



dr 
dx 



v dv 
u c 2 dx 



c 2 



X 



(14) 



However, |fi 2r | 2 x is just the expected Newtonian inertial acceleration, a/v, and dr/dx 
— tan 6* (see figure 6). Hence 



, n a N 

tan 9 — — 



(15) 



An implication of the function r = r(v) is that a full description of the motion of 
bodies in three dimensional space s(x,y,z) necessarily requires consideration of a 5- 
dimensional frame of reference (i.e. s(x, y,z), r and t). The implications of the equation 
tan 9 = ajy/a r are twofold: first, that all finite angular momentun vectors have a 
component in the r dimension, and second that the plane of rotation of a body will 
lie solely in the three space dimensions s(x,y,z) only if the centrifugal acceleration of 
the body tends to zero. One consequence of these principles is that the rest vector 
(defined as Q Q ) for all zero angular velocities must be aligned in the direction of the 5 th 
dimension in the sense of r positive. This is the same as saying that the default zero 
angular velocity component lies in the 5 th dimension for all matter and energy. As the 
total angular velocity, Q, increases the vector becomes increasingly displaced towards 
the ^-direction, i.e. being always normal to the plane of rotation (xyr space in figure 6), 
subtends an angle 9 away from the r-direction. 

Returning to the 'centrifugal' accelerations expressed in a 5-dimensional framework, 
we can now use the expression for tan 9 in equation 1 1 which gives 



— cln 



(a N / a r ) (l + a 2 N / afj +i(l + a%/af) 



-0.5 



(16) 
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which, for small (a^ << a r ), becomes: 

a x = — + ia,N (17) 
a r 

The real part of this expression recovers the MOND limit if a r is identified directly with 
the a Q parameter from MOND. 

For large a at (a^v » a r ), it becomes: 

a x — + i a r (18) 

the real part of which recovers the Newtonian limit. Hence the real terms in these 
equations reproduce the key MOND equations and do indeed argue for a clear association 
between this derivation and MOND. Moreover this analysis provides a form for the 
transition function fi of MOND directly as equation lfjjj]. It can be seen from the above 
derivation that this transition function is governed ultimately by equation 12, which is 
the Lorentz contraction formula approximation for v << c. In addition the value of a r is 
now fixed by equation 5 and in the next subsection we show that its value is remarkably 
close to the value of a Q used by MOND. 

Substitution of equation 18 into equation 3, to satisfy the condition of a force 
balance in 5-dimensions, shows that g x = as expected in the Newtonian limit. 
However if we do the same with equation 17 there is a residual imaginary term which in 
the limit a at — > becomes ia r . This appears to be consistent with a seemless transition at 
the fringes of galaxies from fixed orbits to those increasingly governed by the expansion 
of the Universe. 



4-1. The value of a r 

To calculate a value for a r from equation 5 we adopt a straight-forward cosmology with 
an expanding universe in which time increments, At, and space increments, As can be 
related to shown in figure 7, which gives 

At = -Aa (19) 

c 

As = r u A(psma (20) 

Consider two galaxies embedded in expanding space-time with zero peculiar motion; 

i.e. A0 is constant. Their apparent mutual recession velocity, v r , will be given by ds/dt 

at constant A0 and r u . 

. c . c ctf^ »_ . « 

v r = As — cot a = As — cot — (21) 

1 u 1 u 1 u 

By Hubble's Law this should also be equal to H As which provides the link between H 
and r u 

H= ^ cot ^i (22) 

'V, 

|| This functional form of n, derived here from first principles, is precisely one of the fitting functions 
noted by Milgrom http://www. astro. umd.edu~ssm/mond/faq. html to give good agreement with 
observational data. 
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Figure 7. : Increments of space (As) and time (At) are defined by reference to r. 



If the value for H is taken to be 71 kms~ Mpc~ and the age of the universe, t u , is taken 
to be 13.7 x 10 9 years, then equation 22 gives a minimum value for r u = 6.3 x 10 26 m 
which substituted into equation 5 gives a value of a r — 1.4 x 10 -10 ms -2 . This is in good 
agreement with the MOND parameter a derived from observations of rotation curves of 
galaxies [1,2], so the connection with MOND is now clear and in the remaining sections 
we substitute a Q for aj^jj. If the MOND parameter value of 1.2 x lCr 10 ms -2 is adopted 
then the implied value for r u is 7.5 x 10 26 m. 

4-2. Galaxy dynamics 

The original motivation for MOND, and hence this work, was in explaining two aspects 
of galaxy dynamics; i.e. the flat rotation curves and the Tulley-Fisher relationship. 
Equation 16 allows us to predict what we expect from the derivation in this paper. 
Matching the reals parts of equations 16 and 3 gives 



where M (x) is the mass enclosed within a radius x. The rotation curve, v (x), is then 
given by 



As an example of the transition region for a system with a centralised mass distribution, 
the left-hand panel in figure 8 shows this function for the solar-system. The departure 

This means that equations 4 and 7 can be combined and written as g x = — a ^r, although there may 
be good reasons why a more accurate version is g x — ~a ^- — . 




(23) 




(24) 
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Figure 8. Transition between Newtonian and MOND regimes for the solar system. 



Galactic Orbital Velocities 




Orbital Radius (kpc) 



Figure 9. An example galaxy rotation curve. 



from Newtonian gravity is visually obvious on this plot beyond 300 AU. The asymptotic 
orbital velocity is ~ 360 m/ s. In the right-hand panel we show the difference between the 
Newtonian prediction and that from equation 24. As an example of a galactic rotation 
curve, figure 9 is plotted for a distributed mass inferred from measured data from a 
typical spiral galaxy. In figures 8 & 9 the lower curve show the Newtonian behaviour. 
It is worth noting that in the current derivation there are no free parameters from one 
galaxy to the next, other than needing to know the mass distribution. Uncertainty in 
the value of a Q is a systematic which applies to all galaxies. The asymptotic velocity 
from equation 24 is 

Vtsym = VoGM (OO) (25) 

which underlies the Tulley-Fisher relationship. 



4-3. Linear Motion 

Equation 16 has been derived for a system with angular momentum and predicts an 
observable departure from a Newtonian system. What if we now consider a system 
without angular momentum? Combining equations 9 and 12 we can write 

dr r„ I dll dT\ 
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where U is potential energy and T is kinetic energy of a body with mass m. For a body 

undergoing purely gravitational acceleration towards a much larger mass M without 

any modified response from the Newtonian expectation this becomes 

dr r u ( GM \ 1 

— {-g x ~ a x) = tan0 (27) 



dx c 2 V x 2 ) a Q 
which is none other than equation 1. This implies that there are no components of 
linear momentum vectors lying in the r dimension. Moreover if we also impose an 
energy balance with respect to time, 

into which we substitute equations 9 and 12 we get the result 

9x = a N (29) 

Hence there is a balance condition without reference to the additional dimension 
r (or a a ) and so without departure from Newtonian principles applied in 4-dimensional 
space-time. This implies that MOND-type behaviour is not to be expected from a 
system without angular momentum. 

5. Pioneer Anomaly 

Finally, we consider whether the 5-dimensional orbital dynamics described in this paper 
can be used to explain the Pioneer Anomaly [10]. The acceleration parameter, a r , 
introduced in section 3.1, itself has two components: an angular velocity Q r (with 
magnitude c/r); and c for propagating radiation. The angular velocity, Q r , is orthogonal 
to a r and, being a function of r, will vary as r varies. In particular, Q r will vary with 
time for a body moving away from a gravitating mass, M, in accordance with equation 
9. However, to find the complete time derivative of Q r account must also be taken of 
the time dependency of the radius of curvature of space in an expanding universe as 
described in sections 3.1 and 4.1. Hence we start with the angular velocity Q s that is 
associated with the acceleration parameter, a s , for a space dimension interval As (see 
figures 2 and 10). 

n s = (30) 
r sin a 

Substituting for r from equation 9 (in terms of distance s rather than x) and for a 
from equation 19, the time derivation of Q s can be expressed as 



n 



s 





~r u ~ 




'9s' 


rsma V 


. r _ 




.c 2 . 



ds da\ 

— + cota— 31 
dt dt 



where g s is the gravitational acceleration due to a mass M , ^ is the radial velocity of a 
body with respect to M, (defined as v ra( a a i); ^ is c/r in the general case of r(s) rather 
than r u in equation 19; r u /r ~ 1; and, from equation 22, H = (c/r) cot a, again in the 
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Figure 10. Angular velocities in 5-dimensions 



general case of r(s) rather than r u . The angular acceleration Q r is the component of Cl s 
resolved in the direction orthogonal to a r . Hence equation 31 can be approximated by 

n r = --(g s ^L + Hc) (32) 



where Q r is the angular acceleration derived in the first part (g s v r adiai/c) from motion 
in the direction of s radially away from the M and, in the second part (He), from 
the consequences of unbound motion (i.e. not within a fixed orbit) that is subject to 
the influence of the background expansion of the universe. The angular acceleration 
Cl r is orthogonal to g s and gives rise to a tangential acceleration (i.e an observed 
linear acceleration) acting in the direction of s orthogonal to both Q r and r from the 
relationship 

Tangential acceleration = r A Q r (33) 

This tangential acceleration appears to be Coreolis-like in nature within the context 
of 5-dimensional orbital dynamics - i.e. it derives from angular accelerations that arise 
if there is radial motion (d(r sin a) / dt) orthogonal to an angular velocity (£l s ). It is this 
Coreolis-like acceleration which is identified with the Pioneer Anomaly and is given the 
symbol apA- From equations 32 and 33, apA can be written: 

V radial TT tnA\ 

apA = ~g s He (34) 

c 

Equation 34 can be used to predict values for the Pioneer Anomaly acceleration 
(ap a)- In particular, for Pioneer 11 between distances of 22 AU and 32 AU from the 
sun; and, for Pioneer 10, between 40 AU and 70 AU [11]. The respective average radial 
velocities of the spacecraft at these distances were 11.8 kms" 1 and 12.6 kms -1 . The 
figure for the mass of the sun is taken to be 2 x 10 30 kg and H = 71 kms _1 Mpc _1 . The 
predicted values for apA are shown in table 1. 

A number of things can be noted about these results and equation 34: 
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Distance 
(AU) 


a PA x 10 10 

(ms~ 2 ) 


Pioneer 11 




22 


11.7 


27 


10.1 


32 


9.2 


Pioneer 10 




40 


8.5 


50 


7.9 


60 


7.6 


70 


7.4 


Average 


9.0 



• All the figures for apA lie within the tolerance limits reported for the Pioneer 
Anomaly: 8.74 ± 1.33 x 10~ 10 ms -2 , except for those predicted for Pioneer 11 at 
distances less than 27AU. 

• At the distances considered, the second term (He) makes the larger contribution 
to apA (in fact, 6.9 x 10~ 10 ms~ 2 ) and, in this form of the equation, the term does 
not vary with distance from the sun, nor with radial velocity (v ra diai)- However, 
this term only arises after the Pioneer spacecraft have entered hyperbolic orbits 
(following their final fly-bys) and their motion is subject to the Coreolis-like effect 
described above. 

• If the equation is integrated as a single continuous function over the distance range 
20 AU to 70 AU at an average radial velocity of ~ 12.5 kms -1 and assuming a 
hyperbolic orbit throughout, the value for apA generated is ~ 8.7 x 10~ 10 ms~ 2 . 

• At distances of ~ 15 AU and with the spacecraft travelling at a radial velocity 
~ 10 kms -1 , but without having entered hyperbolic orbit, the value of apA predicted 
by the equation remains ~ 9 x 10~ 10 ms~ 2 as the loss of the second term which falls 
away for a bounded orbit, is made up for by the stronger gravitational acceleration 
of the sun in the first term. 

• As radial velocities decline to zero in bound orbits, so will apA- 

• The equation only takes into account a single space dimension (s) - the radial 
distance from M. This prevents it from describing the full 3-dimensional motion 
of the Pioneer spacecraft as they travel through and out of the solar system. 
In addition, the extent to which the transition from bound to unbound orbits 
progressively introduces the second term is not addressed. A full 5-dimensional 
derivation of the equation should correct this short-coming, but is beyond the scope 
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of this paper. 

• In circumstances where the hyperbolic orbit term does apply, the scale link between 
apA and the MOND parameter ao, is apparent from the relationship He = a cot a 
(derived from equations 5 and 22) where cot a = 5.7 for the current age of the 
universe. 

The level of agreement between the values for apA predicted by equation 34 and those 
reported seems to be reasonably good, albeit that the reported value is interpreted to be 
constant whereas equation 34 predicts a range of possible values that depend on several 
key parameters, such as the radial velocity of the Pioneer spacecrafts. 

6. Discussion 

A model based on simple assumptions about 4-dimensional space-time being closed 
and isotropic, and embedded in a 5 th large-scale dimension r representing the radius of 
curvature of space-time, has been used to derive the balance conditions for a system 
with angular momentum. It has been found that the equations of MOND explaining 
the rotation curves of galaxies appear as a limit within this derivation and that there is 
a universal acceleration parameter, a , with a well defined value, again consistent with 
that used by MOND. In addition this derivation provides, for the first time, a derived 
transition formula between the Newtonian and MOND limits. The following positive 
points should be noted:- 

• Derivation of the Tulley-Fisher relationship from first principles. 

• Derivation of the MOND acceleration parameter (a Q ) from key cosmological 
constants. 

• Identifying Qjq clS du universal acceleration of expansion acting everywhere in 
opposition to gravity (as the universal acceleration of contraction), so that the two 
opposite accelerations jointly define conditions of "tension" within branes occupied 
by matter and energy. 

• Using, as foundation stones, key relationships from General and Special Relativity, 
in particular: 

• Identifying the Lorentz contraction formula as governing the smooth transition 
from MOND to classical Newtonian regimes, as well as from classical Newtonian 
to relativistic regimes - so avoiding the need for an arbitrary MOND transition 
function (//). 

• Generating equations of motion which appear to describe a natural transition from 
fixed orbits to those, on the fringes of galaxies and clusters of galaxies, which are 
affected by the expansion of the universe. 

• A possible derivation of the Pioneer Anomaly using 5-dimensional orbital 
dynamics including a calculation involving key cosmological constants with no free 
parameters. 
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The principal objection to introducing a fifth large-scale dimension into models of 
the universe is generally taken to be a lack of direct observational evidence to support 
its existence. However, the approach adopted in this paper uses a similar premise 
to many other papers in assuming the fifth dimension is quite unlike the three space 
(i.e. visible) dimensions s(x,y,z), in that: first, it has to be treated as mathematically 
imaginary; and second, it does not represent a degree of freedom of motion for matter 
or energy. In consequence, the case for a fifth large-scale dimension must be inferred 
rather than directly observed. Two potential derivations of M OND actually proposed by 
Milgrom - one based on modified inertia and one based on branes in higher dimensional 
space-time - are, to an extent, both reflected in the 5-dimensional model proposed in 
this paper. However, a key conclusion of this paper, which is at odds with these two 
derivations is that departures from classical Newtonian dynamics - as observed within 
four-dimensional space-time - will only be seen in systems with angular velocity or, as in 
the case of the Pioneer Anomaly, where angular accelerations arise within 5-dimensions. 
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